function mom = compute_moments_for_calibration(calib, param, glob, options)

    %% assign parameters

    %    
    param.phi_L   = calib(1); 
    param.sigma   = calib(2); 
    param.epsilon = calib(3)*param.sigma; 
    param.mu_f   = calib(4); 
    param.sigma_f     = calib(5); 
    param.xi      = calib(6);
    param.rho_s = calib(7);

    param.ce       = exp(param.mu_f + param.sigma_f^2/2);


    %% Set up and compute steady state
    
    [param,glob]    = setup(param,glob,options);
    
    eq_SS           = solve_eq_ss(param, glob, options);

    
    %% then compute_moments
    if eq_SS.flag == 0
      mom = compute_moments(eq_SS, param, glob, options);
    else
      mom.m1 = nan;
    end

    mom.flag = eq_SS.flag;
    
end

